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We use current and future simulated data of the growth rate of large scale structure in combination 
with data from supernova, BAO, and CMB surface measurements, in order to put constraints on 
the growth index parameters. We use a recently proposed parameterization of the growth index 
that interpolates between a constant value at high redshifts and a form that accounts for redshift 
dependencies at small redshifts. We also suggest here another exponential parameterization with a 
similar behaviour. The redshift dependent parametrizations provide a sub-percent precision level to 
the numerical growth function, for the full redshift range. Using these redshift parameterizations 
or a constant growth index, we find that current available data from galaxy redshift distortions and 
Lyman-alpha forests is unable to put significant constraints on any of the growth parameters. For 
example both ACDM and flat DGP are allowed by current growth data. We use an MCMC analysis 
' to study constraints from future growth data, and simulate pessimistic and moderate scenarios for 

, the uncertainties. In both scenarios, the redshift parameterizations discussed are able to provide 

significant constraints and rule out models when incorrectly assumed in the analysis. The values 
taken by the constant part of the parameterizations as well as the redshift slopes are all found to 
significantly rule out an incorrect background. We also find that, for our pessimistic scenario, an 
assumed constant growth index over the full redshift range is unable to rule out incorrect models in 
all cases. This is due to the fact that the slope acts as a second discriminator at smaller redshifts 
and therefore provide a significant test to identify the underlying gravity theory. 

Q_|. PACS numbers: 95.36.+x;98.80.Es;04.50.-h 
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I. INTRODUCTION 



Since its discovery over a decade ago [l| , cosmic acceleration stands as one of the most important and challenging 
problems in all physics, see for example the reviews [2| and references therein. As discussed in these reviews and 
others, cosmic acceleration can be caused by the presence of a dark energy component in the universe or alternatively 
\ a modification of gravity physics (namely General Relativity) at cosmological scales. A significant step toward the 
understanding of the cause of cosmic acceleration is to be able to distinguish between the two competing alternatives. 

Indeed, one ongoing approach to understand the origin of cosmic acceleration is to constrain the equation of state 
of dark energy [2J while other approaches rely on comparisons of the cosmic expansion history to the growth rate of 
large scale structure, see the incomplete list |3T-[l5| . Namely, it was shown in many of these references and others that 
two gravitational theories can have very degenerate Hubble curves but yet have distinct functions of the growth rate 
of large scale structure in the universe. 

The growth rate can thus be used in order to constrain the underlying gravity theory and there has been much 
interest in providing a parameterization of the growth factor function with one or two parameters that are distinct 
and characteristic for a given gravity theory, again see the partial list |3Ml5f . 

In this paper, we compare parameterizations of the growth index as a function of the redshift to current and future 
observations. We also introduce a new exponential parameterization for the growth parameter which allows us to 
easily characterize the asymptotic value of the growth index. For current constraints, we use the growth data (mainly 
from galaxy redshift distortions and Lyman-alpha forests) from [l7l - [26j . the Constitution compilation of supernova 
data sets [23], baryon acoustic oscillation (BAO) measurement from the Sloan Digital Sky Survey (SDSS) [29|, and 
the distance to the surface of last scattering of the CMB as measured from the Wilkinson Microwave Anisotropy 
Probe 5 yr data (WMAP5) [31]. We also explore, using a Monte-Carlo-Markov-Chain (MCMC) analysis, how well 
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future growth data from galaxy redshift distortions and Lyman-alpha forests will be able to constrain the growth 
index parameters and to rule out incorrectly assumed underlying gravity theories. 

II. GROWTH OF LINEAR PERTURBATIONS 

The growth rate of large scale structure is derived from matter density perturbation S — Sp rn /p rn in the linear 
regime that satisfies the simple differential equation: 

5 + 2H5-4nG effPm 5 = 0, (1) 

where the effect of modified gravity is introduced via the expression for G e f t . Equation ([1]) can be written in terms 
of the logarithmic growth factor / = din 6/ din a, for many gravity theories as 

/'■ /M~ -^|/ = f%^, ( 2 ) 




where ' is for d/dlna. The growth function / was shown to be well approximated by the ansatz [32| - l35l |: 

/ = n& (3) 

where 7 is the growth index parameter. Reference [32j ] proposed the approximation f(z = 0) = 5^q' 6 for matter 
dominated models and then a more accurate approximation f(z = 0) = ^q^ 7 was proposed in references [33l l34l| . The 
authors of [35] considered dark energy models with slowly varying equation of state and found and index parameter 
7(fi m ,w). This has been the subject of more recent discussions, see for example [H[36j], and also extended to models 
with non zero curvature in [37| and (39J. 

As we discuss below, since growth data spans a wide range of redshift and the growth index is a function of the 
redshift, it is worth exploring parameterization of the growth index as a function of the redshift. For example, in 
references 0, H3| the authors proposed a redshift dependent parameterization of the growth index that was applied for 
redshifts < z < 0.5 and was given as: 7(2) = 70+7' z, where 7' = ^-{z = 0). This showed already the usefulness of 
a variable growth index to distinguish between dark energy models and modified gravity models §,|4(|. Most recently, 
reference [16J proposed a parameterization that covers a wider range of redshift and interpolates to asymptotic values 
where the value becomes almost a constant including redshifts up to at the CMB scale, and we summarize it in the 
next sub-section. 

A. Interpolated redshift parameterization of the growth factor index 

In order to be able to consider constraints from growth data that span over a high redshift range such as for example 
z = — 3.8 growth data provided in [l7l - [26l |. and also future higher redshifts data and constraints up to the CMB 
scale, the authors proposed a parameterization, see (l6| . that interpolates smoothly between low/intermediate and 
high-redshift ranges up to the CMB surface. Similar to the interpolation proposed in |31| for the equation of state of 
dark energy (see appendix C there), the proposed parameterization for the growth index is as follows: 



7( a ) = 7(a) T~TZ 77 + 7„ri„ 1 I — ( 4 ) 



1 1 
h 7 

l + a ttc /a eary l + a/a ttc 

so that 7(a) interpolates between the asymptotic 7 cari!J value at very high redshifts and 

liate ( a ) = 7(a) = 70 + (1 - «ba (5) 

at low redshifts. 

Similarly, using a = 1/(1 + z), the parameterization reads 

= 7(z) : i+ z + 7oo — l +z ^ (6) 

and interpolates between 700 at high redshift (z > zttc) up to the CMB scale and the following form at lower redshifts 
(i.e. for z < z U c) 

l{z) = llateiz) = 70 + 7a (t^Tt) ( 7 ) 
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Parameters for various QCDM models. 


[Wo, W a ) 




7b 


(-0.8,0) 


0.5525 


0.01124 


(-1,0) 


0.5457 


0.01034 


(-1.2,0) 


0.5403 


0.00983 


(-0.8,-0.3) 


0.5460 


0.01512 


(-1.2,0.8) 


0.5579 


-0.00372 


Parameters for some EDE models. 


(w ,C) 


7°o 


7b 


(-0.972, 1.858) 


0.5522 


0.00455 


Parameters for various DGP models. 




Too 


7b 


0.21 


0.6866 


-0.02815 


0.262 


0.6873 


-0.02314 


0.31 


0.6876 


-0.01920 



TABLE I: Best fit values for the exponential parameterization to the growth rate function as generated from the growth 
differential equation corresponding to quintessence dark energy models (QCDM) with constant and variable equations of state, 
early dark energy models (EDE) (see for example [If|), and DGP models. 



We note here that the z Uc is the redshift of transition from a varying growth index parameter 7(2) to an almost 
constant 700 and is not necessarily the same as z t which characterizes the transition from a decelerating cosmic 
expansion to an accelerating one. 

As discussed in [l6j], the growth index for dark energy models (QCDM) using this parameterization takes the form 
given by equation ([6]) where 7(2:) is given by equation ([7]), and where 

QCDM = ^ = WI 

5 — Qw 

As it should, Eq. J5J) reduces to the well-known r y^° M — yy for the ACDM model. 

It was found in [16| that the redshift interpolated parameterization of the growth index provides significant im- 
provements over a constant index value. It was found that the / = QJ n (z) fits the growth function that is integrated 
numerically from the differential equation to better than 0.0065% for models with constant equations of state and to 
better than 0.15% for models with variable equations of state. 

For the Dvali-Gabadadze-Porrati (DGP) models [38], the interpolated parameterization is given by equations (j6]) 
and ([7]) but where the asymptotic value is 0, [36| , 



and as for example derived in Q from 



DGP = 11 fgl 

7oo 16 W 



^iVol^ 1 -^- (10) 



The interpolated parameterization was found in [la ] to provide a fit to better than 0.04% to the growth function for 
DGP model with fit = 0.27 . 



B. Exponential Parameterization for the growth index 

Motivated by the redshift dependence of the growth index, we propose here and explore yet another alternative 
parameterization for the growth index. We use the fact that the growth index should be nearly constant for early 
times and vary only during late times. For this reason we introduce the parameterization for the growth index: 

lexp(z) =700 +7fc(e~ 2/2 «0- (11) 

Looking at the parameterization one can very easily arrive at the relation 7' = —~/b/zttc- This parameterization is 
constructed in such a way one can quickly identify z ttc as the redshift at which the growth index begins to transition 
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FIG. 1: Theoretical fittings for the exponential parameterization. TOP-LEFT: We consider various QCDM models and plot 

the relative error — m -, — - in order to compare the fit of the proposed parameterization to that of the growth factor, /, that 
is numerically integrated from the growth ODE. For ACDM we find the best fit parameters = 0.5457, 76 = 0.0103, 
with z Uc fixed at 0.61. The fit approximates the growth function / to better than 0.015% while the best fit, ^const M — 0.5509 
approximates the growth to 0.65% reflecting the discrepancy at low redshifts. TOP-RIGHT: We plot "fewp{z) = 7oo+7b(e~ z ' /Zttc ) 

for various GR models. BOTTOM-LEFT: We consider the DGP model and plot the relative error "* j. — - in order to compare 
the fit of the proposed parameterization to that of the growth factor f nU m that is numerically integrated from the growth ODE. 
We find the best fit parameters 700 = 0.687, 75 = —0.0231, with fixed z Uc = 0.61. The fit approximates the growth function / 
to better than 0.06% for 0,^ = 0.262 while the best fit, constant j^nst = 0.6795 approximates the growth to 2% reflecting the 
discrepancy at low redshifts. BOTTOM-RIGHT: We plot j e xp(z) = 700 + 7t.(e _Z//z « c ) for various values of f2„, finding for all 
models 700 ~ 0.687, matching theoretical predictions. 



to the constant value 700. It is also worth noting here too, as with the interpolation parameterization above, z Uc is 
not the same redshift, zt, which characterizes the transition from a decelerating cosmic expansion to an accelerating 
one. We must also note that our analysis finds that unless there is an extremely small amount of error on the data 
(on the order of 0.1%) it is difficult to constrain this parameter very well. In other words, there is a large range for 
Zttc for which paramctrization provides an excellent fit to the growth. This is similar to the discussion in WMAP 
paper 31] about the transition redshift for the equation of state where a large range of redshift was found to fit the 
data equally well. For our analysis, we fix zttc to 0.61 which is found to fit a wide range of models very well. 

Fitting parameterization to theoretical data (from numerically integrated growth ODE ([2])) for ACDM and DGP 
models we find that it is able to fit ACDM to within 0.015% and DGP to within 0.09% for values of f2^ raging from 
0.21 to 0.31. It also shows the ability to correctly pick out the theoretical value of 700 for both DGP and ACDM to 
within 0.13% and 0.05%, respectively. The results of the fits are provided in Figure Q] and Table Q] 



5 



2 


fobs 


References 


0.15 


0.49 ±0.1 


[19, 20] 


0.35 


0.7 ±0.18 


[21] 


0.55 


0.75 ±0.18 


[22] 


0.77 


0.91 ±0.36 


(19] 


1.4 


0.9 ±0.24 


[23] 


3.0 


1.46 ±0.29 


[24] 


2.125 - 2.72 


0.74 ±0.24 


[25] 


2.2 - 3 


0.99 ± 1.16 


[26] 


2.4 - 3.2 


1.13 ± 1.07 


[26] 


2.6-3.4 


1.66 ±1.35 


[26] 


2.8-3.6 


1.43 ± 1.34 


[26] 


3-3.8 


1.3 ± 1.5 


[26] 



TABLE II: Summary of observational data used for the growth factor / from galaxy redshift distortions and Lyman-Alpha 
forests as compiled in references [I3-Qj| from original references provided to the right of the table. 



III. OBSERVATIONAL CONSTRAINTS 



A. Current data 



In this section, we will use the parameterizations (j6|) and (fTT|) to fit available observational data considering QCDM 
models and the DGP model. We perform the best fit of the parameters in the various models by % 2 minimization, ie. 



minimizing \ — xtn ^Xbao ~^xi m b y i a a maximum likelihood analysis. In this work, for Type la SNe data, we use the 
397 Type la SNe that compose the Constitution compilation [27[ ■ This collection of samples consists of 250 high-z 
Type la SNe and 57 nearby Type la SNe in the Union set, which were collected by multiple experiments including 
the ESSENCE Survey [4ll E3|. Supernova Legacy Survey p3| . an extended dataset of distant SNe observed with the 
Hubble space telescope, and older observed SNe data, [28| as well as 90 Type la SNe from the CFA3 survey [27j ■ For 
Type la SNe data, we define the x 2 for fitting is as 



X 2 

Asn 



397 

E 

; = 1 



(12) 



where [i(z) = 5 log 10 [di(z)/Mpc] + 25 is the extinction-corrected distance modulus, the total uncertainty in the SNe 
data is Oi. We note that we marginalize over the absolute SNe magnitude, A4, when performing our x 2 minimization. 
The luminosity distance is given by the usual relation: 



(Il(z) 



l + z 



-.S 



dz' 



where 



[ sin( v / |n fe |a;), 



S(y/)Sh\x) = I yW^ 



[sinh(^|fi fe | 



E(z') 



if fife < 0, 
if fife = 0, 
if fife > 0, 



(13) 



(14) 



and the dimensionless Hubble parameter for dark energy with an equation of state parameterized as w{z) — w ± 
w a [z/(l + z)] is given as E(z) = H(z)/H a = {fi (l + zf ± (1 - fi )[(l ± z) 3 ( 1 + w °+ w ») exp(-3w Q z/(l + z))]} 1/2 and 
for the spatially flat DGP model it is: E{z) = (1 - fi )/2 + [fi (l ± zf + (1 - fi ) 2 /4] 1 / 2 . For the ACDM model, we 
use (w ,w a ) = (-1,0). 



To get constraints from observations of the CMB we follow [3l| and define three fitting parameters for comparison 
to WMAP5 data: the shift parameter, i?, 

R(z.) = Vn^H (l + z*)D A (z*), (15) 
with the redshift of the surface of last scattering, z„, given by: 

z* = 1048[l±0.00124(O b /i 2 )-°- 738 ][l±. 9 i(fi m /i 2 ) 92 ]. (16) 
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Constraints on interpolated parameterization for the growth 
index, Eq. J6]), using current data 


Model 


Best fit parameters 






70 






Ho 


X'/dof 


LCDM 


o.92ti:Sg 


— i.4y_ 6 08 


U.Z01_ 023 


72.3±23 


1.19 


DGP 


0.89ii;<f 2 


-4.90 


u.zoo_ 026 




1.29 


Constraints on exponential parameterization for the growth 
index, Eq. (1111). using current data 


Model 


Best fit parameters 






7oo 


7b 




Ho 


x'/dof 


LCDM 


0.16l{;« 


0.59± 2 ;^ 


u.z*y_ 022 


72.51^ 


1.19 


DGP 


o.o3±!i:?? 


0.82ti;- 


n 2t . q +U.U3U 

u.zo»_ 019 


66.O+2? 


1.30 


Constraints on a constant growth index, 7, using current data 


Model 


7 




Ho 


x'/dof 


LCDM 


0.60t££ 






1.19 


DGP 


U - 04 -0.24 


262 +,) ' lm 
u.zuz_ 022 


65.6ii;« 


1.30 



TABLE III: Constraints from the combined current observational data including the growth data given in TABLE [IT] on the 
growth index parameters, SlJ,, and Ho. We see that current data is unable to constrain any of the growth index parameters 
well enough to draw any conclusions. 



The constants 171 and 172 in the above expression are: 

0. 

T+ 39.5(fW) - 763 ' 



0.0783(r>^ 2 )- - 238 

■91 = ; I on c /o LMTia ' ( 17 ) 



and 

l + 21.1(fi 6 /i 2 ) 1 - 81 ' 



0.560 

92 = -, 01 WO 7 9 \ 1 SI ■ ( 18 ) 



A third parameter, the acoustic scale, l a , is 



l. = (l + *)2^W (19) 



r s (z*) 

with the proper angular diameter distance, Da(z) = Dl(z)/(1 + z) 2 and the comoving sound horizon, r s (z*) 

1 da 



r s {z*) = / = , (20) 

v^io a?H(a)y/l + (3r> 6 /4fi 7 )a 

where fi 7 = 2.469 x 10~ 5 /V 2 for T cmb = 2.725K. 

Together the parameters Xi = (R, l a ,z*) are used to fit Xcmb ~ AxiCov^ 1 (xiXj)Axj with Axi = Xi — x° hs and 
Cov~ x (xjXj)is the inverse covariance matrix for the parameters. 

Next, for the BAO, we follow [2!| and define the ratio of the sound horizon, r s (z^) to the effective distance, Dy as 
a fit for SDSS by 

2 f r s (z d )/D v (z = 0.2) -0.198 ^ ra (z d )/JV(z = 0.35) - 0.1094 ^ 

Xbao _ V O0058 J + V O0033 J ' (21) 

with 

Ov(z)- (D 2 A (z)(l + z) 2 j^--y /3 : (22) 



and the redshift, Zd as 

= I • 0.659(O„/i 2 ) -* 



** = , 1 ?!^r/r 9 !n 5 L [i + win (23) 
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0.2 0.3 0.4 0.5 0.6 0.7 0.8 O.J 



FIG. 2: As discussed in section III-B, this is a first simple test to check (using future simulated data) the working of the best 
fit method for the growth index parameters using our MCMC analysis. The growth data is generated using the moderate 
scenario using ACDM model as the background model as well as the fitted theoretical model. TOP LEFT: Interpolated 
parameterization: x 2 /dof — 0.854. TOP RIGHT: Exponential parametrization: \ 2 /dof = 0.855. As expected, for each 
parametrization, we recover best fit parameters of the ACDM model that are practically the same as the ones found from 
other fits in [l(| for the interpolated parametrization, and in Figure [T] of this paper for the exponential parametrization. We 
additionally recover the best fit values for f2„ and Ho that deviate less than 0.4% from the values used to generate the simulated 
future data. 



where 



and 



bx = 0.313(r!„ l ^)- - 419 [l + 0.607(r! m /i 2 ) u - b74 ] 



b 2 = 0.238(ft m /i 2 ) 

We also add the prior Ho — 72 ± 8 km/s/Mpc given by [H 
For the growth factor data, we define 



2\0.223 



(24) 



(25) 



12 

E 



\f obs {z l )-W{z l )^ 
4 



(26) 



where o~fi is the la uncertainty in the f(z) data. 

The growth data that we use (Table [TTJ) is data that was converted in the work of p7l - fl9j from either measurement 
of redshift distorsion parameter j3 — f/b (where b is the bias that measures how galaxies trace the mass density field) 
or from various power spectrum amplitudes from Lyman-a Forest data. The list of the respective original references 
is also given in the table. A caveat in using this data to constrain other cosmological models is that at various steps 
in the process of analysing or converting the data (in the references above and those in Table |TT] ) , the ACDM model 
was assumed. So if one wants to use the data to constrain other models, and in particular modified gravity models, 
one than should redo all the steps assuming that model, starting from original observations. Therefore one should 
keep in mind this caveat when looking at the results below for the DGP model. 

Monte Carlo Markov Chains (MCMC's) are used to compute the likelihoods for the parameters in the model. This 
method randomly chooses values for the above parameters and, based on the % 2 obtained, either accepts or rejects 
the set of parameters via the Metropolis-Hastings algorithm. When a set of parameters is accepted it is added to the 
chain and forms a new starting point for the next step. The process then repeated until the specified convergence 
is reached. In this work the MCMC code we use is based on a modified version of the publicly available package 
CosmoMC [13. 

We perform the fits using ACDM and flat DGP with all the combined data above and allow the following concerned 
parameters to vary: the physical dark matter density, Q,d m h 2 ; the physical baryon density, f^/i 2 ; the ratio of the 
sound horizon to the angular diameter distance, 9 (f2^ and Hq are derived parameters in CosmoMC); and the growth 
parameters. The values found for the growth parameters as well as and Hq are summarized in Table [TTTT We find 
that current data is unable to put significant constraints on any forms of the growth factor index parameters and we 
turn to future simulated data. 
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FIG. 3: Interpolated parameterization. TOP LEFT: Moderate scenario fitting background DGP data to an assumed ACDM 
theoretical model: xV do / = 1-038. TOP RIGHT: Pessimistic scenario fitting background DGP data to an assumed ACDM 
theoretical model: xV rfo / = 1-023. BOTTOM LEFT: Moderate scenario fitting background ACDM data to an assumed DGP 
theoretical model: x 2 /dof = 1.037. BOTTOM RIGHT: Pessimistic scenario fitting background ACDM data to an assumed 
DGP theoretical model: /dof = 1.017. As shown on the figures, in each case the incorrect assumed background model is 
ruled out to 95.45%. 



B. Simulated future data 

The next question that one needs to address is how well future data could constrain the growth index parameters. 
We use here a modified version of the program CosmoMC 47J to include growth rate codes and likelihood function in 
order to find the best fit growth index parameters. We simulate growth data and uncertainties using growth ODEs' 
codes and use in the analysis pessimistic and moderate scenarios based on values of the one sigma uncertainties on 
future growth data f b s - We also simulate future data for Type la SNe, BAO's and CMB surface measurements and 
use all the simulated data sets in order to see what constraints future data can put on the growth parameters even 
when other cosmological parameters are varied. 

To simulate future growth data, we use the discussions and some current data provided in the following references 
[l9l [25l. I2H I49M5II I for data coming from redshift distortions or from Lyman- Alpha forests, in order to extrapolate two 
scenarios. It is expected that the limitations will come from the systematics in each probes. We generate 80 points 
(or bins) for the growth rate that are almost equally spaced by Az = 0.05 between redshifts 1 and 4, but where we 
allow for random small departures from these exact locations up to 1.5a were a is the uncertainty on the simulated 
data point as described below. 

We use a pessimistic scenario where we assume that we will have more data but the uncertainties will get only 
slightly better then the ones of the current data and we use the one sigma uncertainties on the observed values of 
the growth as follows: a(fobs) = 20% for the range < z < 2.0 and cr(/obs) = 30% for 2 < z < 4.0. The moderate 
scenario is: <j{f bs) = 15% (an improvement of a factor of 1.33) for the range < z < 2.0 and a(fobs) = 25% (an 
improvement of a factor of 1.2) for 2 < z < 4.0. 

We simulate the future data by using our best-fit results for cosmological parameters from current data to generate 
the cosmological distances parameters for each future data set, i.e. n(z), / 6 S , r s (zd)/Dv, and (i?, l a , z*) for supernovae, 
growth, BAO, and CMB surface, respectively. For BAO future constraints we used as guidance a Fisher matrix analysis 
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FIG. 4: Exponential parameterization. TOP LEFT: Moderate scenario fitting background DGP data to an assumed ACDM 
theoretical model: xV do / = L 032. TOP RIGHT: Pessimistic scenario fitting background DGP data to an assumed ACDM 
theoretical model: xjdof = 1.022. BOTTOM LEFT: Moderate scenario fitting background ACDM data on an assumed DGP 
theoretical model: x 2 /dof = 1.032. BOTTOM RIGHT: Pessimistic scenario fitting background ACDM data to an assumed 
DGP theoretical model: % 2 /dof = 1 .015. As shown on the figures, in each case the incorrect assumed background model is 
ruled out to 95.45%. 



performed in [52|]. As a conservative approach, we generate seven data points 0.2 < z < 0.7 with a spacing of about 
Az = 0.075 and an uncertainty of 2% for z < 0.5 and 7.5% for 0.5 < z < 0.7 to use for BAO constraints. Here too 
we allow small random departures from the exact data values up to 1.5er. Our future data is in line with predictions 
from BOSS documentation]!!]]. For Type la SNe we generate 500 equally spaced data points (distance moduli) for 
z = 0.0036 — 1.8 with an uncertainty of 0.7%, in-line with uncertainties on these values from future surveys. Again 
we allow small random departures from the exact values up to 1.5o\ We also use the "Bluebook" documentation of 
the Planck satellite|54] to improve the constraints on the CMB surface measurements from WMAP5. Conservatively, 
we decrease the uncertainty by half. 

When we do the fits we use all the simulated data sets and allow for the following parameters to vary: Qdmh 2 , 
flbh 2 , 6, and the growth parameters corresponding to the parameterization used for 7, (as usual, fi^ and Hq are the 
derived parameters in CosmoMC). 

We call the background model in our analysis the model that we use in order to generate the simulated data. 
When we simulate the data for that background model, we use all the functions involved to be those of that model. 
For example, we use the growth equation ([2]) with H(z) and G e // for that model as given explicitly in our previous 
work [l6| for ACDM and DGP models. Similarly, we use the functions of the background model for the supernova 
luminosity distance and for the angular diameter distances for the BAO and distances to the CMB surface. In that 
way all the data has been generated consistently for one background model. 

We call the fitted theoretical model in this analysis the model that we use in order to fit the data previously 
simulated using the background model. In other words, the fitted theoretical model is the assumed theoretical model 
that one will use in order to fit the real data from a real survey. 

Now, one also needs a strategy on how to apply the tests. A trivial exercise is to generate the data using a 
given background model as the true cosmology (ACDM, let's say) and then run the best fit test for j(z) and other 
cosmological parameters using the same theoretical model as the one used to generate the data (ACDM here). In this 
case, one should expect to recover confidence contours where, of course, ACDM should be the best fit (or not too far 
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FIG. 5: Constant Growth Index 7. TOP LEFT: Moderate scenario fitting background DGP data to an assumed ACDM 
theoretical model: x /d°f = 1.063. TOP RIGHT: Pessimistic scenario fitting background DGP data on an assumed ACDM 
theoretical model: % jdaf = 1.041. BOTTOM LEFT: Moderate scenario fitting background ACDM data to an assumed DGP 
theoretical model: X 2 /dof = 1.051. BOTTOM RIGHT: Pessimistic scenario fitting background ACDM data to an assumed 
DGP theoretical model: x 2 /dof = 1.024. 



from it) and the confidence contours sizes will depend on the uncertainties added to the generated data. As shown in 
our Figure [5] this trivial check works very well, recovering best- fit value for growth parameters very close to those we 
would expect. In this case, when these fits are performed we also find best fit values for and Ho that deviate less 
than 0.4% from the values used to generate the data. 

Of course, a more significant strategy that we adopt here is to explore when is the test able to rule out a given 
model. For this, we generate the data using a given background model (ACDM, let's say) and then apply the best 
fit test for 7(0) and other cosmological parameters using a different fitted theoretical model. In this approach, we are 
testing if the simulated data and uncertainties (i.e. the survey) are able to rule out the growth index parameters of 
the fitted theoretical model, as it should, or the uncertainties are too large and more precision is needed from the 
survey. 

In fact, the values of the growth index parameters given by the literature are attained by taking a given model 
and evaluating growth index parameters within that same model and then the characteristic values of for example 
7 = 0.55 for ACDM and 7 = 0.68 for DGP are obtained, see for example [D,[36[. This approach gives values for growth 
index parameters that, by their very construction, are model dependent. For this reason one cannot simply fit growth 
index parameters using an arbitrary model and expect that the values returned will be values indicating what the 
correct model of gravity is. More clearly, for a given background model (i.e. the simulated or true data), the values of 
the growth index parameters that one will obtain from the fit will be different depending on what theoretical model 
is used for the fit. For example, if one is fitting the growth index and assuming the cosmology is ACDM (theoretical 
model used) and finds that 7 = 0.55 is not within the range of uncertainties, it can be concluded ACDM is not the 
correct true cosmology. However, the value measured cannot be compared to the value 7 = 0.68 of the DGP model 
because in this example we are using ACDM as our theoretical cosmology. In other words, it does not matter if 
7 = 0.68 falls within the uncertainties or not because DGP is not the theoretical model used in this example. The 
particular measured value of 7 cannot be used to confirm or rule out a theoretical model other than the one used 
in the fit. For this reason, the growth index parameters are best used to rule out or confirm an assumed theoretical 
model. 

Our results are summarized in Table IIIIB1 Figures [3J H] and [SJ We find that in either the pessimistic or moderate 
scenario for future data uncertainties, the incorrect fitted model is ruled out using either of the two parameterizations 
discussed. We also find that a constant growth index parameter should be able to rule out an incorrect model in the 
moderate scenario but fails in the pessimistic scenario to rule out ACDM when the DGP data simulated data is used. 
In addition to the values of 70 and 700, the slopes of the interpolated or exponential parameterizations (i.e. j a and 
7b) by themselves also rule out an incorrect background. 
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Interpolated parameterization for the growth index, Eq. Q6D. 
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Exponential parameterization for the growth index, Eq. (1111). 
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Constant growth index, 7. 
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TABLE IV: Here we summarize our best fits for multiple parameterizations of the growth index, when we fit generated 
background data to the wrong model as described above. We see that all of the parameterizations are able to find inconsistencies 
in at least one parameter compared to their expected theoretical values given right below. In these fits we find best fit values 
of fi^, ~ 0.230 and Ho — 74.3 when the background model is DGP and the fitted theoretical model is ACDM; and ~ 0.295 
and Ho — 63.4 when the background is ACDM and the fitted theoretical model is DGP. For these values of we would expect 
the following parameter values for the fitted theoretical models: (0.564, —0.025) for ACDM and (0.655, 0.046) for DGP in the 
interpolation parameterization; (0.546, 0.011) for ACDM and (0.688, —0.020) for DGP in the exponential parameterization; and 
(0.552) for ACDM and (0.680) for DGP with a constant 7. 



IV. CONCLUSION 

We explored comparisons of redshift parameterizations of the growth factor index to current and future growth 
data. The first parametrization used was introduced in previous work and interpolates between a redshift dependent 
form at small redshifts and a constant value at high redshifts. A second parametrization based on an exponential 
form is introduced here and exhibits a similar redshift dependence, as it should. We found it to fit theoretical data 
to within 0.015% for ACDM and 0.09% for DGP, over the entire redshift range up to the CMB surface. While 
more precise parametrizations are welcome, we consider that the more significant plus in these redshift dependent 
parametrizations is that they provide the slope of the growth index as a second test to the underlying gravity model. 
This is relevant because the slope is related to variations in j(z) at small redshifts where more data can be obtained. 
Using redshift dependent parameterizations or constant value of the growth index, we find that current growth data 
from redshift distortions and Lyman alpha forests is unable to put significant constraints on the growth parameters. 
In order to explore how well future growth data could constrain these parameters, we simulated growth data and ran 
a Monte-Carlo-Markov-Chain analysis. We find that a pessimistic or moderate scenarios for future data uncertainties 
will be able to rule out an incorrectly assumed theoretical model using any of the two parameterizations discussed 
while we find that in our pessimistic scenario a constant growth index parameter will be unable to rule out an incorrect 
model. This is due to the fact that the slope acts as a second discriminator at smaller redshifts. 
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